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DECUS Progrom Library Write-up 


DECUS NO. 8-556 


CHISQ Chi Square Pro grain 


Ihis program will compute up to 20 chi squares at a time 
on tables as large as 8 by 9. Missing data is permitted. If 
the III is missing or out of range, the entire case is ignored. 

If any item is missing or out of range, that item for that case 
is ignored. 

Input is in the form of identifying information, which 
determines in which row the subject belongs for every chi square 
table; and the individual's scores on up to 20 items, each score 
determining the column in which he beLongs for the appropriate 
chi square table. 'The ID and scores can appear in any order. 

Tables need not all be the same size for the same run. 

Items may have different ranges. 

There is an option for computing the contingency coefficient 
far each table. 
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The program consists of 2 chained segments. The first 
contains a main program and subroutine IRD, and is called 
CHISO* The second is named SQUARE. Each must be compiled, 
loaded, and saved separately on the system device. To allow 
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for execution-time formatting of input data, load CHISQ 
and change the contents of location 4072 from 4124 to 
3400 before saving in the system device. For an 
explanation of execution-time format statements, see the 
accompanying write-up. 
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How to run CHISQ Chi Square Program: 


1) .R CHISQ (Starting Address is 02435) 

'n a heading to appear on output (maximum length is 
36 characters), 

3; Type in number of rows (ID categories). Maximum is 8. 

4) Type in number of items (not including ID). Maximum is 20. 

6) Type “ IOORR ’ contingency coefficient option. (1 is yes , 

0 is no). 

’yie in number of categories for each item (in 2011 format). 
The number of categories can vary from item to item. 

7) Type input device (1-2. 3 or A) t-f •;»■> + • 

v or 4) * If ln P u t is from magnetic 

tape, the file must reside on unit 1 with the name ' CHIN.D1* 

8) Type in which variable is ID. If it is the first, type 1. 
ir it is last and there are 20 items, type 21, etc. 

9) Type in the format statement. 2 lines of 72 characters 
each are permitted. (Be sure to press RETURN twice). 

All data must be read in integer format. The last card 
(or record) must have ID equal to 9. 


Mote: The "HEAD" statement is at 

occur a card reader error, STOP, 
STaRT. 


location 04065. If there should 
load Address 04065, and press 
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Execution-time format Statements 



Execution-time format statements are not permitted by 
the compiler. However, by changing one location in the program 
after it is in core, it is possible to read in a format 


statement (from the teletype only), and use it for reading the 
input. Set up the program (or subprogram) as follows: 

DIMENSION FM'J? (12) 

4 

This is the array where the format statement will 
C be stored. 


15 

C 

c 


read (1,15) FMT 
FORMAT (12A6) 

A format statement as long as 72 characters can 
be read in here. 



20 

C 


25 

C 


READ (IN,25) DATA!,DATA2 


Input data is read in from device IN 
FORMAT (F3.0) 

This is a dummy format that will not be used. 


f 

f 

r 
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t 
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END 
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Be sure to obtain an assembly listing of the program, 
and also a map at execution time, so that the precise location 
to be changed can be determined. 

Load the program and all subroutines, and begin execution* 
When the program pauses or waits to accept data, the appropriate 
location can be changed. 

Now look at the assembly listing. Determine the actual 
address of the first word of array H4T. (See page C-4 of 8K 
&ABR Assembler manual or PS/8 manual for computing actual 
addresses.) Now look through the listing for the Fortran statement 
where the reading of the data is done, in this case, statement 
20. You will notice that in the assembled code will appear 
the following; 


4033 CALI 2, READ 

6201 ARG (n 
DDDD 

6201 ARG /fa 


(call to subroutine READ 
with 2 arguments) 

(n is the device number) 

(DDDD is the address where 
the device number is stored.) 
(fa- is the format number, in 
this case, 25) 


■ (FFFF is the address of the 

beginning of the format 
statement.) 

In our example, FFFF is the address of the fii^st word 
of format statement 25. This is the location to be changed, 
so that it contains the address of the first word of array 
FMT. Therefore, compute the address where FFFF appears and change 
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its contents accordingly. (As a caeck on your computations, 
first compute the actual value 1'FFF, then load the address 
where IfH appears, and see if the contents agree with your 
computation of IWi\ Then reload the address and deposit the 

address for tne first word of FMT.) 

once this one location has been changed, save the core 

image on tape. 


Altering format Statements Already in Core 

format statements are stored 2 characters per word in 
packed, b-oit ASCII code, ihus, the first half of the first 
word of a format statement is always a 50, the stripped ASCil 
code for open parenthesis. F2.G takes 2 locations in storage, 
and appears as 0662 for P2, and 5660 for .0 • ihereiore, one can 
change all or any part of a format statement in core by 
determining its exact location and substituting the codes lor 
the desired format. For any format statement that is likely to 
be changed in this manner, it is advisable to leave plenty of 
extra spaces when writing the program. 
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